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Abstract 

Using the collective field theory approach of large-iV generalized two-dimensional 



Yang-Mills theory on cylinder, it is shown that the classical equation of motion of 
collective field is a generalized Hopf equation. Then, using the Itzykson-Zuber integral 



at the large- N limit, it is found that the classical Young tableau density, which satisfies 
the saddle-point equation and determines the large- N limit of free energy, is the inverse 
of the solution of this generalized Hopf equation, at a certain point. 
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1 Introduction 

The 2D Yang-Mills theory (YM2) is a theoretical tool for understanding one of the most 
important theories of particle physics, i.e., QCD 4 . It is known that the YM 2 theory is a 
solvable model, and in the recent years there have been much efforts to analyze the different 
aspects of this theory. The lattice formulation of YM2 has been known for a long time |3J , and 
many of the physical quantities of this model, e.g. the partition function and the expectation 
values of the Wilson loops, have been calculated in this context [2,3]. The continuum (path 
integral) approach of YM 2 has also been studied in |§] and, using this approaches, besides the 
above mentioned quantities, the Green functions of field strengths have also been calculated 

§• 

It is known that the YM2 theory is defined by the Lagrangian tr(F 2 ) on a Riemann 
surface. In an equivalent formulation, one can use itr(BF) + tr(I? 2 ) as the Lagrangian of 
this model, where B is an auxiliary pseudo-scalar field in the adjoint representation of the 
gauge group. Path integration over the field B leaves an effective Lagrangian of the form 
tr(F 2 ). 

Now the YM2 theory is essentially characterized by two important properties: invariance 
under area-preserving diffeomorphisms and the lack of propagating degrees of freedom. These 
properties are not unique to the itr(BF) + tr(B 2 ) Lagrangian, but rather are shared by a 
wide class of theories, called the generalized 2D Yang-Mills theories (gYM 2 's). These theories 
are defined by replacing the tr(B 2 ) term by an arbitrary class function f(B) ||. Several 
properties of gYM 2 theories have been studied in recent years, for example the partition 
function [7,8], and the Green functions on arbitrary Riemann surface |J. 

One of the important features of YM 2 , and also gYM 2 's, is its behaviour in the case of 
large gauge groups, i.e., the large-iV behaviour of SU(N) (or U(N)) gauge theories. Study 
of the large-iV limit of these theories is motivated on one hand by an attempt to find a 
string representation of QCD in four dimension ||10| . It was shown that the coefficients of 
1/N expansion of the partition function of SU(N) gauge theories are determined by a sum 
over maps from a two-dimensional surface onto the two-dimensional target space. These 



kinds of calculations have been done in |TlJ and [12| for YM 2 and in || for gYM 2 . 

On the other hand, the study of the large-iV limits is useful in exploring more general 
properties of large-iV QCD. To do this, one must calculate, for example, the large-iV be- 
haviour of the free energy of these theories. This is done by replacing the sum over irreducible 
representations of SU(N) (or U(N)), appearing in the expressions of partition function, by 
a path integral over continuous Young tableaus, and calculating the area-dependence of the 
free energy from the saddle-point configuration. In |L3|, the logarithmic behaviour of the 



free energy of U(N) YM 2 on a sphere with area A < A c = n 2 has been obtained, and in [|14[ 
the authors have considered the case A > A c and proved the existence of a third-order phase 

1 



transition in YM 2 . A fact that has been known earlier in the context of lattice formulation 



15| . In the case of gYM 2 models, the same transition has been shown for f(B) = tr(B 4 ) in 



T§ and for f(B) = tr(5 6 ) and f(B) = tr(S 2 ) + gtr(B 4 ) in [0, all on the sphere. 



Such kinds of investigations are much more involved in the cases of surfaces with bound- 
aries. This is because in these cases, the characters of the group elements, which specify the 
boundary conditions, enter in the expressions of the partition functions and this makes the 
saddle-point equations too complicated. In |18| (see also [0J), the authors considered the 



YM 2 theory on cylinder and investigated its large- N behaviour. If we denote the two circles 
forming the boundaries of the cylinder by C\ and C 2 , then the boundary conditions are 
specified by fixed holonomy matrices Uc x =Pexp/ Ci y4 At (x)dx M and Uc 2 =Pexpf C2 A /i (x)dx /i . 
In the large-iV limit, in which the eigenvalues of these matrices become continuous, the 
eigenvalue densities of Uc x and Uc 2 are denoted by o"i(#) and er 2 (0), respectively, where 
9 e [0, 27r]. Then it was shown that the free energy of YM 2 on cylinder, minus some known 
functions, satisfies a Hamilton- Jacobi equation with a Hamiltonian describing a fluid of a 
certain negative pressure |]IB|| . The time coordinate of this system is the area of the cylinder 
between one end and a loop (0 < t < A), and its position coordinate is 9, and there are two 
boundary conditions a(9)\ t=0 = &\{9) and a{9)\ t =A = cr 2 (9). It is found that the classical 
equation of motion of this fluid is the Hopf (or Burgers) equation. Further, it was shown 
that the Young tableau density p c , satisfying the saddle-point equation (and therefore spec- 
ifying the representation which has the dominant contribution in the partition function at 
large-iV), satisfies 7rp c (— 7t<jq(9)) = 9. <Jq{9) is cr(9,t) at a time (area) t at which the fluid is 
at rest. When Uc\ = Uc 2 , 0"o(0) is the solution of Hopf equation at t — A/ 2. In the case of 
a disc, <7 2 (9) = 6(9), the authors have calculated the critical area A c by using the results of 
the Itzykson-Zuber integral jn| at large-iV limit. 



Studying the same problem for gYM 2 has begun in |p0|| , in the context of master field 
formalism. In this paper we study this problem, gYM 2 on cylinder, using the above described 
technique. The plan of the paper is as following. In section 2, by calculating the classical 
Hamilton-Jacobi equation, we obtain the corresponding Hamiltonian for the eigenvalue den- 
sity for almost general f(B) and find the classical equations of motion. It is found that these 
equations are the generalized Hopf equation. In section 3, we show that the Young tableau 
density p c is the inverse function of the solution of the generalized Hopf equation at some 
certain time (area). 



2 The collective field theory and the generalized Hopf 
equation 

As it is shown in [7,8,9], the partition function of a gYM on a cylinder is 



R 



-AC(R) 



where U\ and Ui are Wilson loops corresponding to the boundaries of the cylinder, the 
summation is over all irreducible representations of the gauge group, \R is the trace of the 
representation R of the group element, and C is a certain Casimir of the group, characterizing 
the particular gYM2 theory we are working with. For the gauge group U(N), the group we 
are working with, the representation R is labeled by N integers l\ to l^, satisfying 

h<lj, i< J- (2) 

The group element U (an NxN unitary matrix) has N eigenvalues e ldl to e 10N . The character 
Xr(U) is then 

det(e a 'M 

*« (C7) = -jW' (3) 

where 

j{e ie -} = n(e^-e^). (4) 

j<k 

The Casimir C is a function of /j's. In its simplest form, C has an expression 

N 

C(R) = Y,c(h), (5) 

where c is an arbitrary function. Here, we restrict ourselves to this form. 

In the Large- iV limit, it is convenient to define a set of scaled parameters y^ instead of 
Zi's: 

In the same limit, also two density functions p and a, corresponding to the distribution of 
j/i's and ^j's, respectively, are defined: 

1 N 
°W) ■= nEW-V (7) 

JV 3=1 

P(y) ■= ^E%-%) ( 8 ) 

iV 3=1 



Inserting (H), (f|), and (0) in ([J), using (|6]), and making some obvious redefinition of the 
function c, one arrives at 



R 



i i i L P - NA J2k 9(yk) 

noi}v{el\ 



Here we have defined 



v{9k} ■= n sin 



. Oj — &k 



j<k 



and K is an unimportant constant. 

To proceed, we use the change of variable 



(9) 



(10) 



T k := i9 k 



;n) 



and rewrite fl9|) as 



R 



det [e N ^ } det [e N ^ T t } _ NA y 



D{rl)D{rl} 



g(yk) 



where 



D{r k } := II smh 
j<k 



Tj ~ T k 



We can then use, along the line of [0 , 



Z = Ke 



N 2 F 



(12) 



(13) 



(14) 



and differentiate it to obtain 



OF _ 1 1 ^ 







Ndr k 



(D'Z), 



(15) 



where D 1 is D{r k }. The function g is assumed to be a polynomial. So, to calculate the 
right-hand side of flT5|), let's first calculate it for a monomial. We have 



%n{4A n ^z) 



ND 1 Z^ \Ndr k 



ND'Z £ U 



<9 



iV<9r fc 



D 1 







(n—m) 



Ndn 



1 E" 

fe,m \ / 
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D 1 \Ndr k 



D 1 



:<r+^ 



(16) 



In the large-iV limit, one can of course omit the 0(1/N 2 ) term (which contains higher 
derivatives of F). In this limit, one must also note the limiting behaviours 



k 



dr 



a[r 



)b(l 



(17) 



and 



N 



db_ 

dr k 



d_ 



5b 



Sa(r) 



\T=T k - 



Here we have used a density function <t(t) instead of <j(8), corresponding to the change-of- 



variable (11 



Returning to (|HJ), we define 



D k 



1 d 



In ID 1 ! 



Ndr k 

= y^ coth 

This remains finite, as iV tends to infinity. One then has 



T , 



(19) 




Nlh k 



D 1 




+ D k 



(20) 



Here, the subscript s denotes that part of expression which contains only the derivatives of 
D k , not D k itself. The first equality simply comes from (D 1 )~ l N~ 1 (d/dr k )D 1 = N~ 1 (d/dr k ) + 
D k . To obtain the second equality, one may consider a term with / factors of D k . There are 
m\/[l\(m — /)!] ways to choose / factors of D k from m factors persent. The other D k s, either 
are differentiated or are not present. 

It may seem that this s part vanishes at the large- iV limit, since it contains factors of 
1/N. This is, however, not the case, since there are singular terms in D k at j ~ k. To 
calculate the non-vanishing part of this expression, let us define a generating function: 



QkW 




1 d 



Ndr k 



D k 



(21) 



This is easily seen to be 



Qk{u) 



1 u 9 



D 1 

d 1 (n + j 

D\r k ) 



n 



sinh ( 



Tk-Tj+u/N 



^ k sinh (^) 



(22) 



The s-part of this expression is contained in terms with small values for r k — r,-. To obtain 



this, we use 



Tk - T~ 



k-j 



(23) 



Na(r k y 

let j run from — oo to oo (but j ^ k), and keep only the leading terms in r k — Tj. It 
is easily seen that if one uses this prescription for D k itself, D k vanishes. So, using the 
above-mentioned prescription in (^) gives exactly qk s ( u ), the s-part of q k {u)- That is, 



QksW 



n 

j<k { 



n=l 



k—j , u_ 
Nfffa) ~T JV 



Na(r k ) 

ua(r k ) 




(24) 



So, 



<7fcsW 



sin[7rti5"(rfc)] 



7IW(TfcJ 

Having found this, we return to (p0| ) and arrive at 



*(*£)"*-? (7)" 



r'a,[7ra(r fc )]', 



(25) 



(26) 



where the coefficients a/ are defined through 



SIM 
X 



°' J 



7T 



a; 



z=o "• 
' lcos(vr//2) , 

~ Jb • 



E 

z=o 



Z! l + l 



Inserting ( ^6|) in (|16|), one obtains 

(D'Z) = 



i v n ^ 



(27) 



ND^Z^\Ndr K 



^i:){:y,iMnm--(N d £ 



k,m,l 



h^Mm:~->"'i NS £ 



k.l 



^-t^>'0« + D * 



n—l 



k.l 



2>(?)/*W[*»M]'{!: 



_5S_ 
5a{r) 



n—l 



(28) 



In the last step, we have defined a function S through 



Njf-S := Nj^F + D k 

OT k dr k K 

dS_ d£ 

dA ■ dA 



(29) 



Combining fllS]) with (|28|), we arrive at 
OS 



^^WI^Hr^Hr)]'^ 



5S 

6&{t) 



n—l 



(30) 



where g n 's are the coefficients of the Taylor-series expansion of g. Considering A as a 
time variable, (|30"D can be regarded as the Hamilton-Jacobi equation corresponding to the 
Hamiltonian 



ff = E |«/dra(r)Mr)f 



dr 



-i n—l 



where II is the momentum conjugate to a. 

The summations in ( |3~1~D are easily carried out to yield 



H= 2^ ldTlG 



i-Ka(r) + 






G 



-Z7TO"(r) + 



dU 



where G is an integral of g: 



dG ( \ 
ax 



From (32) , one can obtain the equations of motion for a and II: 



a 



5H 



and 



n 



J_d_ 
2iri dr 



5H 

2m 



( . „, , du\ ( . _ , air 

g I nra{T) + — I - g I -m:a{T) + — 



# «7tct(t) 



Defining 



v := 



dr) ' 

dU 



g -nra[T) 



0lT 
7fr 



(31) 



(32) 



(33) 



(34) 



(35) 



(36) 



as a velocity field, in correspondence with what defined in [18,19], one can combine fl34|) and 
(|35| ) into a generalized Hopf equation: 

d 
v ± iixo = — — — [g(v ± 27Tcr)] (37) 

In the case of YM 2 , where g(y.k) = \y\, this equation reduces to Hopf equation found in 



|18| . In this case, when one of the boundaries shrinks, so that one has a disc instead of a 
cylinder, that is o~ 2 (#) = S(6), the Itzykson-Zuber integral can be used to obtain a solution 
for the Hopf equation and from that the critical area of the disc has been obtained [I8|. In 
our problem, gYM 2 , we do not know such an integral representation. 

3 The dominant representation 

It is shown in jnj that the character Sr([7), can be written as 

Xr(U) = e N2s ^\ (38) 

where, for large N, 

~ [p, or] = S[p, a] + \ J dy y 2 p(y) + B[a]. (39) 

Here B is some functional of a, and £ satisfies a Hamilton- Jacobi equation 



lr = 5 /d ''" ( -'' ) 



' a sx\ 2 i 2 2 , 



(40) 



in which the variable /x satisfies the boundary conditions 

li(x,t = ti) = a(x), (41) 

and 

/i(x,t = l) = p(x), (42) 

Defining V as the derivative of the momentum conjugate to /i, as in the previous section, 
it is seen that the Hamilton-Jacobi equation ( |40|) is equivalent to the following evolution 
equation for /i and V. 

a-*&=°- (43 > 

where 

<& :=V + inn. (44) 

and V = |-^. 

Using (|39| ) in (|1]) and (|I2|), we see that in the large- N limit, the dominant representation 
satisfies the following saddle-point equation 

8 



^hm = ^- (45) 

or 

where the summation is over boundaries. Now recall fl3"T|). Rewriting it as 

f~+-^9(f)=0, (47) 

where 

/ : = v - ma, (48) 

one can write an implicit solution to (|47"D as 



f(r,b)^f{r-(b-a)g'[f(r,b)],a}, (49) 

where a and b are two particular values of the time variable (here, actually the area). The 
same thing can be done to solve (fl3|). In fact, one can define two functions H + and H~ as 

H + (x):=x-(t-T)$(x,T), (50) 

and 

H-(x):=x+(t-T)®(x,t), (51) 

and see that $ is a solution to fl43|) if these two functions are inverses of each other, i.e. 
H-[H+(x)} =x p|. 



As an ansatz for H ± (with t = 1 and T = 0), we take 

H^ := A#' + iffp, (52) 

and 

#+ := fl* - oto-j 

= /* (53) 

where 

f(x)=f {x-A i g'[fi(x)}}. (54) 

Here A* is the area between a curve for which / is /o and the i-th boundary. The meaning 
of (|5^) and fl53|) is that we are seeking a solution to (|3) with boundary conditions 



Vi(x, 1) = Ai5r'(x) - x, (55) 



and 

/i(x, 0) = a(x). (56) 

For such a solution, (flop is obviously satisfied. 
It is now easily seen that 

Hr[Ht{x)\ = A t g'[f t (x)] + inp(f {x - A ig '[f t (x)]}). (57) 

If fo(x) is the inverse of iirp, 

inplf (x)} = x, (58) 

then we have 

Hr[H?(x)]=x, (59) 

So PS) is satisfied if p satisfies (|5lf). But / is generally a complex function, whereas p 



should be real. It is now better to return to the earlier variables 9 and a. We have 

ia(i6) = a(6). (60) 



So (^f) can be written as 

np[-ifo(x)) = x, (61) 

where 

/ := v — ina. (62) 

The argument of p in (|6T| ) is real if t>o is zero. So, if there exits a loop on the surface, for 
which the velocity field is zero, then there exists a dominant representation p, satisfying 

7rp[-7RT (#)] = 6. (63) 

This is the same as that obtained in [I8| . Note, however, that the equation governing the 



evolution of a, (p^), is different from the corresponding equation in [|T^ . 
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